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I. Introduction 


Wireless receivers process signals that bear information as well as disturbances caused by 
the transmitter/receiver circuits and channel impairments such as fad:ing, interference, and 
additive white gaussian noise (AWGN). Usually, the receiver knows only some 
statistical properties of the signal and disturbances. From these statistical properties and 
using a finite observation of the received signal, the receiver is able to estimate the 
transmitted data symbols. The receiver makes the decision on the received data using 
locally generated symbol clock and carrier oscillator, both of which are not referenced to 
the actual versions used to generate the data at the transmitter. The receiver has to 
estimate the offset between locally generated carrier and symbol clock to those used at 
the transmitter. Clock mismatches are labeled as time jitters, while local carrier 
mismatches are labeled either as phase rotation errors or frequency offset errors. Carrier 
phase recovery is the process known as carrier tracking, where the receiver estimates the 
offset between the local oscillator phase and the actual phase of the transmitted carrier. 
Carrier frequency offset recovery is known as carrier acquisition, and it is the process of 
estimating the offset between the frequency drift/change of the local oscillator and the 
actual (received) carrier frequency transmitted at the transmitter. The methods in which 
the receiver estimates carrier mismatches can be feed forward (open Loops) or feedback 
(closed loops) synchronization. Feed forward method is characterized by the absence of 
hang up and cycles slips, and fast and accurate acquisition. Such advantages are absent 
with conventional feedback synchronizers such as Costas loop. Another advantage of 
feed forward algorithms is that the synchronizer can be fully implemented using digital 
signal processing techniques as shown in figure 1. 
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Figure 1: Complete Feed Forward Baseband synchronization system 


This paper reports an overview of existing as well as presenting new carrier acquisition 
and tracking systems that are used in modem satellite and personal communication 
receivers. In particular algorithms that acquire and track the carrier frequency in the 
absence of symbol timing or a data preamble sequence. The objective is to investigate 
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these algorithms. Then verify their applicability and performance in realistic situations 
such as wireless hand held cellular and satellite TDMA receivers. The performance 
evaluation of the algorithms is conducted through simulations and theoretical analysis. 
The work is performed in two parts. In the first part; we report the mathematical 
development of these algorithms and their statistical properties and “pull in” range in 
AWGN channels. In the second part, we present the performance analysis of the 
recovery process through computer simulation for a variety of applications concerning 
satellite communications receivers, in particular, techniques proposed that acquires 
carrier offsets when the initial error (offset) is too large to employ conventional loops 
(open loops or closed loops). Furthermore, we present a practical means of implementing 
these algorithms using commercial digital signal processor chips. 


II. Carrier Frequency Acquisition and Tracking 


In addition to removing carrier frequency offsets, phase errors between the received 
signal and the local oscillator has to be removed as well. Phase error can result in 
rotation of the received symbols resulting in degradation of the bit error rate and 
ultimately degradation in the receiver sensitivity. The BER deterioration as a result of 
phase mismatches is defined as being the amount of increase in signal power (in dB) that 
is needed in order to get the same BER performance from a receiver that uses an ideal 
carrier synchronizer [ 1,2,3,4,5,6,7,8,9], and given by (1) 


D = 4.3(1+ E, /N&T; [dB] (1) 


where E, /N, is the symbol energy to noise ratio, o° is the phase variance of the 
estimated phase out of the synchronizer. For second order Costas loop, the phase variance 
is given by (2) 
% = s Fad’ 1 (2) 

where B, is the one sided loop bandwidth, and 7 is the symbol rate (loop iteration rate). 
Acquiring and tracking carriers using Costas loops for short-burst satellite TDMA signals 
can prove to be a difficult task. The loop acquisition time as a function of loop bandwidth 
and the frequency offset to be acquired [6] is given by (3) 


y 
_n (2xAf) 27 > “Reel (3) 
acq 16<(2~B,)~ 2nrB, 

Wherex< is the second order type loop damping factor and Af is the frequency offset. To 
achieve faster acquisition, the loop bandwidth in equation 3 has to be increased, while 
doing so, the phase error variance increases in equation 2, leading to reduction of BER in 
equation 1. The term on the right-hand side of equation 3 is the tracking time [6], while 
the term on the left-hand side is the “pull-in” or acquisition time. To demonstrate the loop 
bandwidth, acquisition timing and BER degradation, let us assume a frequency offset due 


to Doppler as 1 kHz and a < = 0.707, figure 2 shows a plot of t,, , phase variance and 
BER degradation as a function of loop bandwidth B,. Clearly as the loop becomes 


narrower to achieve lower BER degradation by having lower ol, the acquisition time 
becomes impractical for typical short TDMA satellite bursts [S]. Meanwhile, if the loop 
bandwidth is widened to achieve faster acquisition, both phase variance and BER 
deterioration become large. Ideally and based on our lab experience, a typical Costas 
loop bandwidth of 25-50 Hz is required for satellite modems with d.ata rates of less than 
56 kbps using MPSK. This allows for a maximum phase variance of 10° rad° (or a 
standard deviation of gz =I.8 1OS°), which results in a BER degradation less than 0.1 dB. 
One way to speed up feedback loops is by utilizing frequency sweep or variable loop 
bandwidth, both of which are complex to utilize and may result in destabilizing the loop 
operation. Loops with a variable loop bandwidth have been utilized [ 121, however their 
implementation creates a multiple magnitude of implementation problems, and needs 
continued maintenance. DFT aided loops have also been used, however a large number 
of DFT bins (high resolution) is required to estimate the carrier to an offset less than 10 
Hz [ 81, such resolution requires huge computational resource and have rendered useless 
for mobile low power receivers. 
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Figure 2: Loop bandwidth versus carrier phase variance and BER deterioration 
for Feedback type carrier recovery systems. 
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11.1 Open Loop M"™ Power Frequency Acquisition Algorithm 


To derive the mathematical representation of this algorithm, we start by representing the 
received matched filter complex M-PSK baseband signal in AWGN n(t) by 


2(t) = ef OOOH) 4 n(p) (4) 


where 4(t) is the phase due to the transmitted symbols, e(t) is the phase term that 
includes the carrier phase error and phase rotation introduced by the propagation channel, 
which may vary in time, however for most of our analysis its assumed to vary slow when 
compared to the variation in phase due to symbols, Af is the frequency offset introduced 
by either channel effects, or offset errors due to the various up and down conversions that 
the signal encounters while being transmitted and received. To perform frequency offset 
detection, we pass the matched filter signal through a differential digital detector which 
eliminates the carrier phase error rotation, by z(kT)z* ((k - N)T) , this yields a signal that 
contains frequency error offsets given by (5) 


2kT)z* ((k - N)T) = expG(2ntkT 
- 2niq((k - N)T) + B(KkT) 

- B((k - N)T) + @(kT) 

- 4((k - NY?) + n’(kW 


where k is the index of the baseband samples, n’(t) is a complex AWGN with real and 
imaginary parts that are independent and have a variance of a,’ = [2E, / NoJ-’ , or 


4 =T[T,2E,IN,]-‘,  whereT, is the sample time. The differential phase detection 


causes cancellation of B(t) , because such phase variation is assumed to be equal between 


i(B(kT)-B((k-N)T)) . : 
two consecutive symbols, so the component e! OES ND m (5) yields a unity term. 


For simplicity we will lump the frequency offset terms due to offset difference between N 


i(2nwT-2nAf((k-N)T)) = ej 
consecutive samples as being é (2nwl2nAf(NIT)) = ej(2WW (notice that we assume 


N=I or one symbol period unless otherwise indicated, or the over-sampling factor, such 


as denoting every 4th sample is a symbol period). The term e@°*''°°**'' °°’ IS due to 
symbol modulations, this term can be removed by processing the signal through an M- 
power non-linearity, where the value A4 is an integer equal to the modulation symmetry 
angle (1.e., for QPSK M=4, for BPSK M=2). QPSK modulation is removed due to the 
4j(i7fn) 

fact that e = 1, forn=0,1,.2,3.... Doing so yields a noisy sine-wave signal in 
equation 5 with a fundamental frequency that is equal to the frequency offset associated 
with the received baseband signal. For QPSK, where M=4, this is given by (6) 


. = e4j(2GtT) ra n’(kT) (6) 


Designating the Ke complex sample of the differential baseband signal as being a real part 
Re(z, ) = I, and an imaginary part Im(z, ) = Q, , making the current matched filter 


sample Z, = I, + jQ, , and the previous sample, which is N-samples back in time is 
Z*M = I+ - je,-, . Substituting both terms into (5) yields (7) 


* 


zizke ji = LET yey 4 Q&-N]+ .&dlr-,v - Qk-iv& 1 (7) 


R&;-N 3 Im(li, ) 


Equation 7 represents a frequency-offset detector. It can be shown that 
LA + jQk-NQk = cos(27~MAfl) and QKIKVN - jQ,-, , = sin(2xMAfl) 


Both terms are commonly referred to as being the dot and cross product terms, 
respectively. Expanding both terms further yields an initial frequency acquisition value 
given by (8) 


2nMAfl= arg {(zkzi-N)M (8) 


where arg{x) = tan-’ (x). In order to reduce random variations of (8) due to AWGN and 
phase residuals from data modulations, equation 8 is averaged over L-symbols yielding a 
final carrier offset estimator algorithm represented by (9) 

Yinl,23.0)"] [Hz] (9) 


= s tan”! k=l 


on yRe[(z,z7.y)"] 
k=l 


This algorithm has a maximum capture range of Af 5 + R, /2M Hz, or roughly 10% of 
the transmitted data rate, where R, =1/7. For instance, using a 10 ksps symbol rate 
satellite modem, and a QPSK modulated signal (A4=4), the upper limit of this algorithm 
is Af= +2500 Hz as shown in figure 3. The algorithm’s estimate accuracy as a function 


of variable SNR is shown in figure 4 for detecting a 1 kHz offset. Clearly the algorithm’s 
performance is vTy accurate for SNR greater than 10 ~—s GB. 
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Figure 4: M" power carrier offset estimator performance for an offset of 1kHz 


After frequency offset correction, the matched filter signal still has an arbitrary phase 
error that have to be estimated and compensated. Before we introduce the phase 


estimation stage, it is of interest to present several other frequencies offset algorithms that 
differ from (9) in term of | implementation and performance. 
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Figure 5: Viterbi carrier offset estimator algorithm 


11.2. Viterbi Frequency Acquisition Algorithm 


Another form of the M” power algorithm was developed by Viterbi [6]. This algorithm 
employes two separate non-linearity operators, one operates on the argument 


A4 arg(&-hl s, and one operates on the magnitude of the Matched filter output 


é 
(z,2%-»|) - This algorithm have been proven to yield better performance for some 


particular modulations for low signal to noise ratios (especially for QAM modulations 
[6]). The Viterbi frequency offset estimates is defined by (10) 


a R fo * "  jM \arg(z,2; 
af = A ast elon (10) 


where d, is a smoothing filter designed to provide robustness when operating at low 


SNR and to counteract slow signal variations that may be present due to multi-path 
attenuation. Figure 5 shows an implementation of this algorithrn. The smoothing filter 
can be a simple single low pass filter with a digital transfer function given by 


diy es ~, where p is a constant such that 0 < p < 1. 
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II.3. Symbol Correlation Based Acquisition Algorithm 


Although simpler to implement, the M" power and Viterbi algorithms presented above 
suffer from self-noise generated by the use of the non-linearity process for removing data 
modulations. For QPSK, the use of the 4" order non-linearity results in a 6-dB rise of 
noise floor, and 3 dB of noise increase for the BPSK signals. Several other algorithms 
have been employed which are based on the sample to sample correlation function and do 
not use non-linearity processing. These algorithms involve more computational 
complexity. The complexity comes from the fact that it takes N(2L—-N-1)/2 complex 
multiplications to perform the auto-correlation of two consecutive complex samples of 
the matched filter. One such algorithm was developed by Luise and Reggiannini in [14]. 
The frequency-offset estimations is given by (11) 


Af = ay" wkJare R(kK)-argR(k-H] (LL) 
T =I 


where the window function is given by 


3[(L -kXL—-k +1)-N(L-N))| 
Wh) = aN? ONL GE? -1 


and the symbol correlation function given by 
L 


1S . 

R(k)= +See (k-J) 
Td kaj 

wherel < j<N, N<L/2. The algorithm's capture range is given by Af = 2R, /10 Hz, 

which is comparable to the Viterbi and the M* power algorithms. 


Kay introduced a similar version to this algorithm, where the frequency estimates is given 
by (12) 
ri N 
gpk’ ¥ w(kJarg{e(k)z" (k -} (12) 


27 ‘=I 


with the window function given by 


ede 6k 


N(N +1)QN +)) 
Luise&W [10] introduced another form of this algorithm, where the frequency estimate is 
given by (13) 


21 


Where the window function (parabolic ) is given by 


aA N~l . « 
Af = R, el mtbermierrat} (13) 
k=0 


6k(L —k) 
k)= TE =1) 
L&W's capture range is Af = R,/10 Hz. 


Fitz introduced another version in [10], with frequency estimation is given by (14) 


et ee 
f= inayat (14) 


with a capture range of Af = R,/2N Hz. 


Finally, L&R introduced another version that is similar to (12), where the frequency 
estimate is given by (15) 


Af = ae yanel SRO) (15) 


with a capture range of Af = R,/N [Hz]. 


Mengali and Morelli reported simulation [10], which showed that these algorithms have 
the same performance for SNR more than 14 dB for QPSK modulations. However, with 
different capture ranges as indicated above. 


11.4. DFT-Based Acquisition Algorithm 


Another technique for estimating the frequency offset is based on DFT [9]. Unlike the 
algorithms discussed before; this algorithm does not require symbol-timing estimates. 
Figure 6 shows the digital implementation of this algorithm. The frequency-offset 
estimate is by (16) 


Af = : anf SIC, (oc, (6 [Hz] (16) 


ut k-0 


.~ 27m 


where C_,(k) = > 2k" (k - Nye’ v ,and C,,(k) = > z(kk' (k- Nye * 


L L 


Notice that C_,(k) and C,,(k) are conjugate pairs, this helps in saving computational 


: R 
resources on the DSP chip. This algorithm has a maximum capture range of Af < i Hz, 


where N is the number of samples per symbols. We have found [10] that this algorithm 
does not yield acceptable performance, unless the matched filter signal has been filtered 
using a bandpass filter. The filter bandwidth varies depending on the roll off factor (@ ) 
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used. The filter should be designed with center frequency at the spectral component of 
interest in z(kT), which contribute to the periodic component of interest in 


z ((k-N)T )z(kT ). Roughly, the BPF should have a center frequency of R,/2 Hz, and 
a bandwidth that is equal to @R,. 
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Figure 6: Joint symbol timing and carrier offset estimation algorithm using DFT 


III. Conclusion 


feedback and feed forward techniques were presented for carrier frequency acquisition. 
Tradeoffs between these algorithms were introduced in term of their implementation 
complexity, their estimation accuracy, and their frequency capture range. Computer 
simulations were used to show their performance in AWGN channels for QPSK 


modulations. 
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